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ABSTRACT 

We investigate the spatial clustering of galaxies in the PSCz galaxy redshift survey, 
as revealed by the two-point correlation function, the luminosity mark correlations, 
and the moments of counts-in-cells. We construct volume-limited subsamples at dif- 
ferent depths, and search for a luminosity dependence of the clustering pattern. We 
find no statistically significant effect in either the two-point correlation function or the 
mark correlations and so we take each subsample (of different characteristic luminos- 
ity) as representing the same statistical process. We then carry out a counts-in-cells 
analysis of the volume- limited subsamples, including a rigorous error calculation based 
on the recent theory of Szapudi, Colombi and Bernardeau. In this way, we derive the 
best estimates to date of the skewness and kurtosis of IRAS galaxies in redshift space. 
Our results agree well with previous measurements in both the parent angular cata- 
logue, and in the derived redshift surveys. This is in contrast with smaller, optically 
selected surveys, were there is a discrepancy between the redshift space and projected 
measurements. Predictions from cold dark matter theory, obtained using the recent 
semi-analytical model of galaxy formation of Benson et al, provide an excellent de- 
scription of our clustering data. 
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1 INTRODUCTION 

The clustering pattern of galaxies observed today reflects 
an interplay between two fundamental processes: the grav- 
itational growth of primordial density fluctuations and the 
physics of galaxy formation. In general, we expect the pro- 
cess of galaxy formation to segregate the galaxies from the 
underlying dark matter distribution and even to give rise 
to dependencies of clustering on physical properties, such as 
galaxy colour, luminosity, morphological type or star forma- 
tion rate (e.g., Kauffmann et al. 1999, Benson et al. 2000). 
Such segregation is usually referred to as "biasing." 

Biasing was originally introduced as a useful device to 
match theoretical predictions for the distribution of cold 
dark matter (CDM) with observations of the spatial correla- 
tions of clusters and galaxies (Kaiser 1984, Davis et al. 1985, 
Bardeen et al. 1986). Since then, the statistical properties of 
the dark matter distribution in CDM models has been ex- 
tensively investigated by means of N-body simulations (e.g. 
Jenkins et al. 1998, Gross et al. 1999). Observationally, how- 



ever, little is known about the large-scale distribution of dark 
matter, although this situation may soon change as gravita- 
tional lensing techniques become increasingly sensitive (Fis- 
cher et al. 2000). In the meantime, partial information on 
biasing may be obtained by studying the relative clustering 
of galaxies with different physical characteristics. This, in 
turn, requires large surveys like the one we analyze in this 
paper. 

The simplest statistical tool to quantify clustering is the 
two-point correlation function, defined as the excess prob- 
ability above (or below) random that an object be found 
at a certain separation from another, randomly chosen ob- 
ject. Two-point statistics completely describe a Gaussian 
point process. However, the galaxy distribution is patently 
non-Gaussian, as evidenced, for example, by the presence of 
rich clusters and walls. A natural generalization is the com- 
plete set of A'^-point correlation functions which provide a 
full characterization of a distribution. Unfortunately, their 
measurement and interpretation become exponentially dif- 
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ficult as the order, N, increases. This is why counts-in-cells 
analysis, which extracts the average of the iV-point corre- 
lation functions over a coll, has become the most practi- 
cal and popular tool for describing higher order clustering 
in the galaxy distribution (e.g., Peebles 1980, Efstathiou et 
al. 1991, Bouchet et al. 1993, Szapudi, Meiksin, & Nichol 
1996, Gaztanaga 1994, Szapudi & Szalay 1997, Szapudi, & 
Gaztaiiaga 1998). 

This paper calculates two-point statistics and presents a 
counts-in-cells analysis of the PSCz survey of IRAS galaxies 
(Saunders et al. 2000). Our analysis is motivated, in part, by 
a desire to characterize some of the possible manifestations 
of biasing, such as a dependency of clustering on luminos- 
ity or a particular functional form for high order statistics. 
To this aim, we construct volume-limited catalogues from 
the PSCz survey which provide a controlled, homogeneous 
framework for statistical analysis. The most efficient tool 
to search for a luminosity dopeudoucc of clustering is the 
luminosity mark correlation (Beisbart and Kerscher 2000) 
which we apply to the PSCz. We also compare our measure- 
ments with a model of galaxy formation in the context of the 
CDM cosmology (Benson et al. 2000) which makes specific 
predictions for the statistics that we investigate. A previous 
search for luminosity effects in the clustering of IRAS galax- 
ies, based on the QDOT survey, failed to reveal any signal 
(Moore et al. 1994). However, a preliminary analysis of the 
PSCz survey (Maddox et al. 2000) suggests that some effect 
could well be present. 

The rest of this paper paper is organized as follows. In 
§2, we present brief details of the dataset and describe our 
measurements of the two-point correlation function, mark 
correlations, and the moments of counts-in-cells up to fourth 
order. In §3, we compare our estimates with measurements 
for other surveys, both optical and infrared, and with the 
prediction of the Benson et al. model. 



2 DATASET AND MEASUREMENTS 

In this section, we characterize the clustering of galaxies 
in the PSCz redshift survey using three complementary de- 
scriptors: the two-point correlation function (§2.2), the mark 
correlations (§2.3), and the higher order correlation ampli- 
tudes or cumulants, the skewness and kurtosis (§2.4). All 
three statistics were calculated for a series of volume-limited 
subsamples, each containing all PSCz galaxies in spheres 
of radius 50, 75, 100, 125, 150, 175, and 200ft-^Mpc, namely 
1160, 2159, 2189, 1985, 1671, 1456, and 1259 galaxies, respec- 
tively. We chose to work with volume-limited subsamples 
because these are homogeneous catalogues that are simple 
to treat statistically and, as shown by Colombi, Szapudi and 
Szalay (1998), the full series of catalogues yields results that 
are equivalent to a measurement with optimal weights in the 
full catalogue. 

2.1 The PSCz Catalogue 

PSCz is a redshift survey of all the IRAS galaxies in the 
Point Source Catalog that have a flux at 60 ixm greater 
than 0.6 Jy. With its 15,500 galaxies distributed over 84.1% 
of the sky, it constitutes the deepest and densest all-sky red- 
shift survey to date. The median depth is just 8100 km s~^, 
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Figure 1. The two-point correlation function of the PSCz survey 
in redshift space. Three-sided symbols represent measurements 
using the Landy-Szalay estimator, while four-sided symbols show 
estimates derived from scaled counts-in-cells (c.f. §2.4). Results 
for three representative volume-limited subsamples arc displayed. 
Statistical crrorbars arc typically smaller than the symbols (ex- 
cept on the largest scales), and arc not plotted for clarity. The 
solid line is the fit by Fisher et al. (1994, see also Seaborne et al. 
1999). The dotted line shows the predictions of the CDM model 
of Benson et al. (2000) for spiral galaxies. 

although useful information is available out to 30,000 km 
s~^ at high galactic latitudes and to 15,000 km s^^ else- 
where. Full details of the survey may be found in Saunders 
et al. (2000). The analysis was restricted to the region of sky 
covered by the survey, with no attempt to interpolate the 
galaxy distribution within the masked areas (Branchini et 
al. 1999). Since the sample is relatively local, we can ignore 
evolutionary effects on the clustering and biasing. 

2.2 Correlation Function 

We estimate the two-point correlation function of the PSCz 
survey by means of the Landy and Szalay (1993; LS) es- 
timator which has been shown to be optimal (Szapudi & 
Szalay 1998; Kerscher, Szapudi, & Szalay 1999). Denoting 
the number of data-data, data-random, and random-random 
pairs in a bin at distance r by DD, DR, and RR respectively, 
the definition of the estimator is 

^{r) = {DD - 2DR + RR)/RR. (1) 

For volume-limited surveys, uniform weighting yields mini- 
mum variance. 

The three-sided symbols in Fig. 1 show our estimate for 
the PSCz two-point correlation function using the Landy- 
Szalay estimator on volume- limited samples of radii 50, 100, 
and 150/1"^^ Mpc. The square symbols show a different esti- 
mate based on the counts-in-cells (CIC) analysis described in 
§2.4. While a rigorous error calculation for the LS estimator 
is not available at present, errors for the CIC estimator can 
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be computed using the FORCE package (see §2.4). The 1-cr 
uncertainties are typically less than 5 — 7%, but grow rapidly 
at large separations, becoming 20 — 25% at r ~ 25/i~^Mpc 
and larger beyond this. The LS and CIC estimators are only 
expected to agree in the regime where the correlation func- 
tion is well approximated by a power-law regime. In this 
regime the uncertainties in both methods are most likely 
comparable. (The error bars have been omitted from the 
figure for clarity.) It is important to recognize that the data 
points in Fig. 1 (and in most other plots in this paper) are 
not independent and so the estimated errors cannot be used 
for maximum likelihood estimation without further analysis. 
This requires computing the cross-correlation matrix of all 
measurements, a complicated and ill-understood procedure 
which is beyond the scope of this letter. 

Results for the remaining volume-limited subsamples 
are consistent with those displayed in the figure. For com- 
parison, the solid line shows the fit of Fisher et al. (1994) to 
the correlation function of the IRAS 1.2-Jy survey, ^(r) = 
(r/ro)^, with ro = 4.43/i"^Mpc, and 7 = 1.28. Our results 
are consistent with this fit, perhaps with a cut-off on large 
scales, as well as with the fit of Seaborne et al. (1999) to 
the full (flux-limited) PSCz survey. Since the deeper sam- 
ples contain predominantly brighter galaxies, the similarity 
of the correlation functions in the subsamples suggests that 
any luminosity dependence of clustering is small. However, 
the different samples have galaxies in common and so their 
correlation functions are not independent. This weakens the 
sensitivity of this test for a luminosity dependence and so 
we prefer, instead, to use the mark correlations, as we now 
describe. 



2.3 Luminosity Mark Correlation Function 

The elegant tool of mark correlations, introduced into astro- 
physics by Beisbart and Kerscher (2000), is ideal for quanti- 
fying the luminosity dependence of clustering. If rrii denotes 
a mark, or a parameter characterizing a class of objects, the 
joint probability of finding a pair of galaxies at (mi, ri) and 
(m2,r2), respectively is 

T{m\,m2,r\,r2)dm\dm2Cpr\d?r2 (2) 

(see also Peebles 1980). If the mark is discrete, dm\ is 
a Stieltjes-Lebesgue measure; otherwise it is the usual 
Lebesgue measure. The two-point correlation function is ob- 
tained by integration over the marginal distribution. 



(1 + = / Tdmidm 



(3) 



where r — \r-i — r2\- The conditional probability of finding a 
pair with marks mi and m2 respectively is 



P(mi, m2|r) 



r(mi,m2,ri,r2) 

n2(i + ew) ■ 



(4) 



This quantity may be estimated from the catalogue as the 
ratio of all the pairs with the prescribed marks divided by 
all the pairs, both numbers taken at a given separation. 

The luminosities were divided into luminosity quartiles, 
each of them containing equal numbers of galaxies. Thus, 
marks m = 0, 1, 2, 3 represent bins of increasing absolute lu- 
minosity. The rest of this subsection is exclusively concerned 




4*m, 




Figure 2. The conditional probability P(mi, m2\r). This is plot- 
ted as a function of 4mi -|- m2 (with m\ > m2) in order to sep- 
arate the points. Here, rrii = 0,1,2,3 represents the luminos- 
ity quartiles of a volume-limited PSCz subsample with radius of 
SO/i^'^Mpc, with the corresponding luminosity ranges between 
(9.0610**2, 1.1823310*3, 1.64892 lO^^^, 2.7935 10*3, 4.2884 lO**) 
ergs/sec. For clarity, only symmetrical errorbars, estimated 
from a series of 10 realistic mock PSCz catalogues, are 
shown. Results for several scales are plotted, with small hor- 
izontal shifts to the right (for clarity), in increasing or- 
der: 2.35, 3.26,4.51,6.25, 8.67, 12. Ol/i-^Mpc (upper panel) and 
16.64, 23.06, 31.94, 44.26, 61.32, 84.95h-iMpc (lower panel). 



with the shallowest volume- limited subsample, cut at a dis- 
tance of 50/i~^Mpc, since this contains the largest range of 
luminosities. The same analysis was repeated for the other 
volume-limited subsamples with similar results. 

The probability, P(mi,m2|r-), is shown in Fig 2, as a 
function of 4mi -I- m2, mi > m2 for a number of pair separa- 
tions. For a particular pair (mi,m2), results corresponding 
to increasing scales are slightly shifted towards the right for 
clarity and only symmetric errorbars are plotted. These were 
calculated from the dispersion of a set of 10 mock PSCz cata- 
logues constructed from an N-body simulation by Cole et al. 
(1998) of a flat ACDM model, with cosmological constant, 
A = 0.7, and without any luminosity segregation. Fluxes 
were assigned randomly to galaxies to match the PSCz lu- 
minosity function and the various selection criteria discussed 
by Branchini et al. (1999) were applied in order to obtain re- 
alistic mock catalogues. For the smaller distances considered 
(upper panel), the errors agree with a naive Poissonian esti- 
mate, while for the larger scales (lower panel), the dispersion 
in the mock catalogues signiflcantly exceeds the Poissonian 
expectation. Therefore, on small scales the dominant source 
of error is discreteness or shot noise, while on large scales 
finite volume effects are important. 

The dotted line in Fig. 2 corresponds to the theoreti- 
cal expectation for the probability, p = (1/4)^ = 0.0625, for 
no correlation between the bins. Points above or below this 
line indicate a positive or negative correlation respectively. 
The results in Fig. 2 are consistent with no correlation on 



© 0000 RAS, MNRAS 000, 000-000 



4 I.Szapudi 



any scale, for all pairs of marks. To quantify this statement, 
we use the 10 mock catalogues just mentioned. The formal 

for all the points plotted, ~ 0.381664, would be un- 
realistically small for 120 independent points. However, the 
points in the figure are not independent because the nor- 
malization condition for P{mi,m2\r) requires the sum over 
all the bins to add up to unity. For this reason, the formal 

has a non-standard distribution which can, in principle, 
be derived by Monte-Carlo methods using our mock cata- 
logues. Although 10 mock catalogues are not sufficient to 
determine the distribution in detail, they allow us to reject 
a significant signal in Fig. 2: the for the sample is close to 
the median for the simulations, 0.381914, whose highest and 
lowest 10 values are 0.204304 and 0.512983, respectively. 

Our results agree with those from a similar, but inde- 
pendent, analysis of the PSCz survey by Kerscher and Beis- 
bart (private communication). However, they disagree with 
preliminary results by Maddox et al. (in preparation) who, 
using different subsamples and statistical techniques, find a 
small but significant dependence of clustering on luminosity 
in the PSCz catalogue. Nevertheless, on the basis of analysis 
performed here, we may assume that any luminosity effects 
are sufficiently small that different volume-limited subsam- 
ples can be combined together to obtain minimum variance 
estimates of the higher order moments, as we do in the next 
section. 

Before proceeding, we note that a comparative study of 
the IRAS 1.2 Jy and PSCz surveys by Teodoro et al. (1999) 
shows that the two density fields differ, within ~ 80/i~^Mpc, 
by a monopole term, a finding which is consistent with the 
spherical harmonics analysis of the PSCz survey by Tadros 
et al. (1999). Both these studies show that this apparent 
discrepancy disappears at a flux limit of 0.75 Jy. We have 
therefore repeated our analysis using samples limited at 0.75 
Jy. Our results remain essentially unchanged for this sample, 
with a formal ~x' = 0.468862 for the comparison in Fig. 2. 

2.4 Counts-in-Cells Analysis 

We carried out a counts-in-cells analysis following the pre- 
scription developed by Szapudi et al. (1999). In a nutshell, 
we measured counts-in-cells using the successive convolu- 
tion algorithm described in that paper, with a high rate of 
oversampling (up to 10^ cells) on all scales for each volume- 
limited subsample. In practice, only about half of the cells 
intersected the geometry of the catalogue embedded in a 
cubic grid. From the counts-in-cells, the skcwncss and kur- 
tosis were calculated using the method of Szapudi & Szalay 
(1993). This technique automatically corrects for shot noise 
by replacing moments with factorial moments. Thus, con- 
tinuous definitions of the cumulants suffice: 



Sn = N^'-'^Qn 



in 



(5) 



where S denotes the field of density fluctuations, and ^ = 
{5^) {N = 3,4 for this paper). In the power-law approxi- 
mation, the average correlation function is proportional to 
the correlation function. A Monte Carlo integration of the 
power-law flt to ^(r) given in the previous subsection yields 
^/2.4 = 4. We plot 4(r) derived from the flrst of equa- 
tions (5) and this relation in Figure 1 (fllled squares). The 
result is consistent with the direct estimate of ^(r) at pair 





Figure 3. Near optimal composite measurements of S3 (upper 
panel) and 54 (lower panel). Different points can come from dif- 
ferent volume-limited cuts, as detailed in Table 1. The errorbars 
were calculated using the SS model. The solid lines with the en- 
veloping dotted lines representing the variance arc the predictions 
of the semi-analytic model of galaxy formation of Benson et al. 
(1999) with the corresponding variance. 



separations at which the power-law approximation is valid. 
For larger scales the above relation is no longer valid. 

Wc determined the S3 and 5*4 for all the volume-limited 
subsamples and, for each cell-size, we cherry-picked the value 
with the smallest error. The analysis of the preceding sec- 
tions suggests that any luminosity bias is small, so we can 
assume that the different volume-limited subsamples rep- 
resent the same statistical process. Colombi et al. (1998) 
have shown that in these circumstances, selecting the high- 
est signal-to-noise measurements from a series of volume- 
limited surveys is equivalent to the optimal minimum vari- 
ance weighting of the flux-limited survey. 

The errors were determined ab initio using the FORCE 
(FORtran for Cosmic Errors) package based on the the- 
ory of Szapudi & Colombi (1996), and Szapudi, Colombi & 
Bcrnardcau (1999). For each volume- limited subsample, the 
error was obtained for each order and scale under a range 
of plausible assumptions for the parameters involved. When- 
ever possible, these parameters were measured from the sur- 
vey itself, such as its volume, the average number-counts and 
the variance in a cell. The variance over the survey volume 
was derived assuming the standard CDM power spectrum 
with r = 0.5 and as = 1. The higher order cumulants for 
3 < S3 < 8 were calculated from perturbation theory with 
the effective index of the power spectrum n = —0.5. This 
yields slightly higher Sns than are typically measured in 
the PSCz survey, thus providing a conservative overestima- 
tion of the (statistical) errors. Finally, the assumption of 
Szapudi & Szalay (1993, SS) was employed for the structure 
of cumulant correlators, the two point analogs of the Qjvs, 
namely Qnm = QnQm (see SS for details). 
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Table 1. Estimates of S:i in cubic cells of size I. The errors, 
""Sn 1 were calculated using the FORCE package. The depth of 
the volume-limited subsample picked out by the near optimal 
procedure (see text) is denoted by R and given in h-^Mpc. The 
value of R corresponding to each estimate of S3 is also applicable 
to the corresponding estimate of ^4, except for the case indicated 
by an *, for which R = 50?i~^Mpc for 54. 
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Provided that the parameters arc well tuned, the calcu- 
lations typically provide better than 50% accuracy for the 
statistical errors, normally a conservative overestimate as 
verified by Colombi et al. (1999) and Hoyle et al. (1999) 
from N-body simulations. It should, however, be borne in 
mind, that if the plausible guesses made for the parame- 
ters and models are grossly incorrect, or if systematic errors 
dominate, the theoretical error estimates might be inaccu- 
rate. We have varied some of the assumptions and param- 
eter values within reasonable bounds without finding sig- 
nificant change. When the relative errors arc approaching 
unity, the perturbative approach breaks down, but still in- 
dicates that the measurement has low significance. For more 
details of the method and its applicability see Szapudi et al. 
(1999). Note that these error estimates do not take into ac- 
coount cross-correlations or possible systematic errors and 
could therefore underestimate the true uncertainties. 



3 DISCUSSION AND SUMMARY 

We have characterized the clustering properties of galaxies 
in the PSCz survey of IRAS galaxies by means of the two- 
point correlation function, mark correlations, and the mo- 
ments of counts- in-cells. Our two estimates of the correlation 
function, one based on a direct measurement and the other 
derived from counts-in-cells, arc consistent with each other 
and agree well with previous measurements for IRAS galax- 
ies (Moore et al al 1994, Fisher et al. 1994, Seaborne et al. 
1999). Neither the correlation function itself nor the mark 
correlations show any significant evidence for a dependence 
of clustering on luminosity, over the limited range of lumi- 
nosities probed by our volume-limited subsamples. Finally, 
we have firmly established that the skewness of the distri- 
bution of counts-in-cells in redshift space for IRAS galaxies 
has a value S'3 ~ 2. 



Our estimated values of Sn, for N = 3,4, are listed 
Table 1 and displayed in Fig. 3. The variation of these 
quantities over the range of scales probed by our analysis 
is small. Over the entire range, S3 = 1.89 ± 0.62 and S4 = 
7.00 ± 4.13, while over the restricted range 1 — 20/i.~^Mpc, 
S3 = 1.93 ± 0.35. These results are in good agreement with 
previous analyses of counts-in-cells of IRAS galaxies, based 
on the QDOT survey (Saunders et al. 1993). For the parent 
catalogue of the 1.2Jy survey, Meiksin et al. (1992) obtained 
S3 = 2.2 ± 0.2 and 54 = 10 ± 3, while for the 1.2 Jy red- 
shift survey itself Bouchct et al. (1993) found S3 = 1.5 ±0.5 
and 54 = 4.4 ± 3.7 and Fry & Gaztafiaga (1994) derived 
S3 = 2.1 ± 0.3 and S4 = 7.5 ± 2.1. The significance of our 
study lies in the more densely sampled data set, the use of 
state-of-the-art measurement techniques, and a rigorous er- 
ror calculation. All these features have enabled us to extend 
the dynamic range of previous studies and achieve unprece- 
dented accuracy. 

In contrast to the situation at optical wavelengths (see 
Hoyle et al. 1999 for an up-to-date discussion), the skewness 
of IRAS galaxies in redshift space agrees well with that in the 
parent, projected catalogue, as measured by Meiksin et al. 
(1992). Presently, these are the only datasets large enough 
to allow a measurement of S3 with a precision better than 
10% over a large dynamic range. The 2dF and SDSS surveys 
will enable estimates of S3 with an error of only a few per- 
cent and of S4 with an error of 10-25% over a large dynamic 
range (Szapudi et al. 1999). These surveys will clarify the 
role of redshift space and projection effects in the apparent 
disagreement between angular and 3D cumulants for opti- 
cally selected galaxies, and enable a more direct comparison 
between the statistics of IRAS and optical galaxies. 

In general, the flatness of the S3 curve in redshift space 
is well understood from theoretical arguments, and sup- 
ported by A''-body simulations. Although the skewness rises 
to a non-linear plateau in real space, the random veloci- 
ties associated with the 'finger of god' distortion in red- 
shift space act as an effective smoothing. As a result, the 
5*3 curve remains flat in redshift space. Naturally, biasing 
complicates this picture inferred from simple theory and 
dark matter simulations. It is therefore perhaps surprising 
that S3 should be so similar for IRAS and optically-selected 
samples since these have somewhat different spatial distri- 
butions. We speculate that this agreement is the result of a 
cancellation effect whereby the differences in spatial distri- 
bution are compensated for by differences in the strength of 
the redshift space distortions acting on each galaxy type. 

To address the effects of biasing on statistical measure- 
ments of the galaxy distributions requires a detailed the- 
ory of galaxy formation. The solid line in Fig. 3. shows the 
predictions of the semi-analytic model of galaxy formation 
proposed by Cole et al. (2000) and Benson et al. (1999). Ac- 
cording to the model, the higher order moments in redshift 
space for spiral galaxies are very similar to those for the 
galaxy population as a whole, although they differ in real 
space (Baugh, Szapudi, & Benson 2000). The predictions in 
the figure refer to spiral galaxies which, as a class, arc a rea- 
sonable representation of IRAS galaxies. The dotted lines 
give the variance computed by varying the input parame- 
ters in the error estimation procedure. The agreement with 
the PSCz results is remarkably good. The shape of the S3 
and S4 curves in the model reflects the kind of cancellation 
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effects mentioned above (see Baugh et al. 2000 and Hoyle 
et al. 1999 for further discussion). As shown in Fig. 1, the 
model predictions also agree well with the two-point corre- 
lation function of the PSCz survey. This model of galaxy 
formation assumes a flat, A-dominated CDM cosmology in 
which galaxies form by hierarchical clustering. The agree- 
ment with our results is the more remarkable since there are 
no adjustable parameters in the comparison, so that the the- 
oretical lines in Figs. 1 and 3 are to be regarded as genuine 
predictions of the model. 



Acknowledgments 

We thank the rest of the PSCz team for use of the data 
before publication, S. Cole for providing mock PSCz cata- 
logues and A. Benson and C. Baugh for providing theoretical 
predictions from the Durham semi-analytic model of galaxy 
formation and for useful comments. IS was supported by the 
PPARC rolling grant for Extragalactic Astronomy and Cos- 
mology at Durham. CSF acknowledges a Leverhulme Re- 
search Fellowship. The FORCE (FORtran for Cosmic Er- 
ror s) package is available from its authors, S. Colombi and 



IS (http:/ /www. cita.utoronto.ca/~szapudi/istvan. html) 



REFERENCES 

Benson, A., Baugh, C, Cole. S.. Frcnk. C .S. & Lacey, C. 1999, 



MNRAS, submitted { istro-ph/9910488) 



Baugh, CM., Szapudi, I., Benson, A. 2000, in preparation 
Beisbar t. C. K erscher. M 2000, ApJ, submitted (astro- 

ph/ jooiosel ) 

Bernardeau, F. 1994, ApJ, 433, 1 
Branchini E., et al. 1999, MNRAS, 308, 1 
Bouchet, F.R., et al. 1993, ApJ, 417, 36 

Cole, S., Hatton, S., Weinberg, D. & Frenk, C.S. 1998, MNRAS, 
300, 945 

Cole, S., Lacey, C, Baugh, C. & Frenk, C.S. 1999, MNRAS, 
submitted. 

Col ombi. S.. Szapudi I.. Szalay, A.S., 1998, MNRAS, accepted 



(astro-ph/9711087) 



Efstathiou, G., et al. 1990, MNRA S. 247. 10 

Fischer et al. 1999 ApJ, submitted( |astro-ph/9912119 ) 



Fisher, K.B., Davis, M., Strauss, M., Yahil, A., Huchra, J., 1994 

MNRAS, 266, 50 
Fry, J. N. & Gaztanaga, E. 1994, ApJ, 425, 1 



Hoyle, F.. Szap udi. I., & Baugh, CM. 1999, submitted (astro 



ph/S 911351) 



Juszkiewicz, R., Bouchet, F. R., & Colombi, S. 1993, ApJ, 412, 
L9 



Kersche r, M., S zapudi. I. & Szalay, A. 1999, submitted (astro- 
ph/ j912088| ) 



Landy, S.D.,& Szalay, A. 1993, ApJ, 412, 64 
Meiksin, A., Szapudi, I., & Szalay, A., 1992, ApJ, 394, 87 
Moore, B., Frenk, C.S., Efstathiou, G. & Saunders, W., 1994, 
MNRAS, 269, 742 



Saunders et al. 2000, MNRAS, (submitted astro-ph/000111 7) 



Seaborne et al 1999, MNRAS, (accepted astro-ph/9905182 J 
Szapudi, I., & Colombi, S. 1996, ApJ, 470, 131 
Szapudi, I., Colombi, S., Bernardeau, F., 1999, MNRAS, accepted 
Szapudi, I. & Szalay, A. 1993, ApJ, 408, 43 (SS) 
Szapudi, I. & Szalay, A. 1998, ApJ, 494, 41L 
Tadros, H. et al. 1999 MNRAS, 305, 527 



Teodoro, L. et al 1999 MNRAS, submitted, (astro-ph/9908358) 



